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We examine a two-dimensional electron waveguide with a cut-circle cavity and conical leads. By 
considering Wigner delay times and the Landauer-Biittiker conductance for this system, we probe 
the effects of the closed billiard energy spectrum on scattering properties in the limit of weakly 
coupled leads. We investigate how lead placement and cavity shape affect these conductance and 
time delay spectra of the waveguide. 
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fQ . I. INTRODUCTION 

O : 

i The scattering properties of ballistic electron waveguides have been studied by a number of authors, but primarily 
in the semi-classical regime, where resonances strongly overlap 0, Q. Several experiments have shown that the 

, shape of the waveguide cavity can strongly affect the conductance, depending on whether the shape induces integrable 
or chaotic motion In this paper, we investigate the scattering properties of an electron waveguide whose scattering 
cavity strongly resembles that of a quantum dot. We investigate this system in a regime where the resonances are well 
^-j ■ separated in energy, and we look at the effect of distortions in shape and changes of lead placement, which change 

, the symmetry of the scattering cavity, on the behavior of the scattering resonances. 

' T ne geometry of the waveguide we consider in this paper is shown in Fig. [I]. The scattering system consists of a 
cavity with two cone shaped leads attached, which we will treat as infinitely long. The walls of both the cavity and 
the leads are assumed to be infinitely hard (infinite potential) walls. The cavity has the shape of a circle of radius 
R, but with a segment of the wall, subtended by an angle u>, replaced by a flat segment which we call the cut. The 
opening angles of the leads are A6*i and A82 ■ The angular positions of the center of the cut and the center of the 
lead openings are denoted by C, and #2, respectively. 

The effect of conical leads has already been studied in earlier works (e.g., Ree and Reichl ||, Persson et al. |J, and 
Berggren, Ji, and Lundberg Qj). Compared to straight leads (which have most commonly been used in publications 
on waveguides), conical leads are more similar to the electron sea which provides the source of electrons in most 
experiments. As shown in Ref. j^], conical leads allow tunneling resonances through the waveguide cavity for energy 
regimes in which conduction is prohibited for straight leads. 

A classical particle moving in a closed billiard whose shape is that of a circle with a cut (the shape of our cavity) 
will display a rich range of chaotic behavior as the size of the cut is varied [|| . Quantum signatures of the classical 
chaos in this billiard have been studied in Ref. ||. One question we will be interested in is how well the scattering 
process can probe this dynamics of the closed cavity. In the low-energy regime we will consider, we expect that the 
waveguide will exhibit well-separated Fano resonances at energies close to the energy eigenvalues of the closed billiard. 

Beyond simply allowing to probe eigenenergies of the closed billiard as scattering resonances, lead placement itself 
can alter the symmetry properties of the closed billiard. We will study the effect of symmetries by investigating how 
the positions of the leads affect the behavior of the Fano resonances. 

In Sec. 2, we will discuss cavity states and the scattering states that are used to obtain the scattering S-matrix for 
a system with conical leads. In Sec. 3, quantum properties of the closed circular cavity with a cut will be reviewed. 
In Sees. 4 and 5, we will look at the behavior of conductance and Wigner delay times for waveguides with full-circle 
and cut-circle cavities. Leads will be attached to cavities symmetrically or asymmetrically. Finally, in Sec. 6, we 
will summarize our results. A short discussion of the numerical method used for our computations is given in the 
Appendix. 
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II. SCATTERING WAVEFUNCTION AND S-MATRIX 



Matter waves inside the waveguide pictured in Fig. [T] are governed by the Schrodinger equation, which, in polar 
coordinates, takes the form 
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Here, k — y2mE/fi 2 is the wavevector of the particle wave, m is the particle mass, Ti is Planck's constant, and E 
is the energy. The geometry of the system, which determines its scattering properties, is taken into account by the 
boundary condition that the scattering wavefunction '5 vanish on the walls. 

The solution 9 is computed numerically by expanding it into suitable basis functions, and determining the expansion 
coefficients from the Schrodinger equation and its boundary conditions by the method descri bed i n the Appendix. 
Below, we will discuss the expressions for 9 in the different regions of our system. In Sec. [I A , we look at the 



wavefunction inside the cavity. The wavefunction inside the leads is considered in Sec. II B, where wc will also define 
the scattering S-matrix of the waveguide. 



A. Cavity states 

In the interior region of the cavity, the solution of Eq. ([!]) can be written as 

9= C ^A- (2) 

7=— n 7 + l 

Here, the / 7 's are suitable expansion functions. The c 7 's are the expansion coefficients, which must be determined 
from the boundary conditions, by the procedure described in the Appendix. The cutoff number ro 7 has to be chosen 
in a proper way to achieve optimal numerical accuracy within a reasonable computation time. We used n 7 = 90 for 
all computations presented here. 

For the functions / 7 , several choices are possible in principle. For our geometry, it proved to be numerically most 
efficient to use solutions of the Hclmholtz equation [(V 2 + k 2 )f 1 = 0, without boundary conditions] which can be 
separated into polar coordinates, {p,f). These functions take the form 

/>,¥>) = J hl (kp)e^, (3) 

where Jui denotes the Bessel function of order |7|. 



B. Lead states 



For waveguide scattering problems, it is usual to express the properties of the system in terms of an S-matrix, which 
contains reflection and transmission amplitudes from incoming channels to outgoing channels in the leads. Before we 
can properly define the S-matrix for our system, we first have to specify what we mean by incoming and outgoing 
channels (also called "modes" ) in this particular case. 

Following the conventions from scattering problems with straight leads, we require that an incoming (outgoing) 
mode x~ (x + ) is a solution of 

(V 2 + k 2 ) X ± =0 (4) 

which is separable into coordinates longitudinal and transversal to the lead walls. As a further condition, the transverse 
part has to be a standing wave vanishing on the walls and the longitudinal part describes a wave propagating towards 
(away from) the cavity. 

For the conical shape of our leads, the longitudinal and transversal coordinates are just the polar coordinates p and 
ip, respectively. Thus, the separation ansatz with a standing wave in the transverse direction (a — 1, 2, . . .) takes the 
form 

x£ = Nitf(p) sin 

where 6^ locates the angular position of the center of the zth lead, R < p < oo, and 6^ — < (p < ftj + A6>j/2 

(i =1, 2 for leads 1 and 2, respectively). A is a (yet unknown) constant factor. 

Inserting this into Eq. (0) leads to an equation for the radial part of the wavefunction, 
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This is the Bessel differential equation and therefore u l ^(p) can be expressed as a linear combination of Bessel 
functions, Jan/A8i(kp) an d Y^/^g^kp). Since we required that the longitudinal part represent a propagating wave, 
we write the radial solution in terms of Hankel functions, 



ut (P) = (kp) = (kp) ± iY^ (kp) 



(7) 



which approach exponential functions at infinity: 
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Here, we used the notation if+ for the Hankel functions of first (second) kind, of order 7. 

The normalization constant, |iV|, is calculated from the condition that each channel carry unit current, i.e., we 
require that 
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Here, the longitudinal component of the probability current density, j l± = e p + e v , is given by 
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This yields 
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With these results, the expression for the incoming/outgoing channel a = 1, 2, . . ., corresponding to the ath standing 
wave in transverse direction, becomes 
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To fix the overall phase factor, N/\N\, in Eq. (12), we require that the solution be real at the interface between 
the leads and the cavity. Then, if a mode is completely reflected back at this position, the corresponding reflection 
coefficient will be 1. This can be achieved by multiplying the modes with a complex phase factor. At p = R, the 
Hankel functions (the only complex factor in x^) can be written as 



iPL (kR) = (kR) 
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Because we want the probability amplitude at this position to be real, we will use the complex phase factor 
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Then, the final form of the incoming and outgoing propagating modes in channel a becomes 
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As we now have derived the expression for incoming and outgoing channels, we are finally in the position to define 
the S-matrix of our system. The matter wave \& inside the leads is expanded in terms of the incoming and outgoing 
modes, x^iPi'-P) an d xL + (/ 9 ' < ^)> respectively. Then the expansion coefficients describe the probability amplitudes of 
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reflection or transmission of the matter wave from an incoming channel to the outgoing channels. More precisely, if 
the electron is incident in lead 1 and channel 0, the wavefunction in the leads takes the form 

*/s,i =xp~ + J2r aax l a + and = Yl ta ^ + ■ ( 16 ) 

a a 

Here, ^0,1 is the state in lead 1, ^0,2 is the state in lead 2, and r a g {t a p) is the probability amplitude that the 
electron enters in channel j3 and is reflected (transmitted) into channel a. Likewise, if the electron is incident in lead 
2 and channel 0, the wavefunction in the leads takes the form 

Vp,2 = xfr +J2r' af3 xl + and ^ = ^ t af}X l + ■ (17) 

a a 

Here ^/3,2 is the state in lead 2, is the state in lead 1, and r aB and t a3 are the reflection and transmission 

probability amplitudes, respectively. Note that the solution '5 of Eq. ([!]) is not unique, but we obtain different 
solutions \P p depending on which channel (3 is assumed to carry the incoming part of the matter wave. 

The S-matrix connects the amplitudes of the incoming propagating electrons to the amplitudes of the outgoing 
propagating electrons and it is therefore constructed out of all the reflection and transmission amplitudes. Thus, the 
S-matrix can be written 

S =(?R')< (18) 

where R and R' are square matrices composed of reflection amplitudes, r a p and r' a g, and T and T' are square matrices 
composed of transmission amplitudes, t a p and t' a g. The reflection and transmission amplitudes must be computed 
numerically (see Appendix). 



III. CLOSED CUT-CIRCLE BILLIARD 



One aspect of the scattering properties of the waveguide that we wish to explore is how close a correspondence 
there is between the Fano scattering resonances and the eigenstates of a closed version of the waveguide cavity. In this 
section, we therefore describe some features of these eigenstates as the cut size on the cut-circle billiard is changed. 

The Schrodinger equation [Eq. (|Tj)] for a quantum particle in a closed full-circle billiard is separable, with one part 
describing the radial and the other part the azimuthal motion. The (discrete) energy spectrum can be labeled with two 
good quantum numbers. One of them, n = 1, 2, is associated with radial motion, and the other, I = 0, ±1, ±2, 
with the angular motion. The cigenfunctions are 

*I,n = ^1 (%|,» P) e llv , (19) 

where Jm is the Bessel function of order \l\ and (k\i\ jTl R) is its nth zero (i.e., Jm vanishes at p = R). The corresponding 
eigenenergies E\i\ n are then calculated via E = h 2 k 2 /(2m). The energy eigenstates of the circle billiard with angular 
momentum |Z|>1, corresponding to clockwise or counter clockwise rotation of the matter wave, produce two- fold 
degenerate eigenenergies. Only the eigenenergies with I = are non-degenerate. 

If we apply a cut to the closed system, the classical dynamics of the particle exhibits hard chaos as long 
as w<180°. Looking at the quantum properties, we observe a splitting of the degenerate eigenenergies as the cut 
is inserted. The absence of degenerate eigenvalues (level repulsion) is a well-known quantum feature of classically 



chaotic systems (see, e.g., Ref. |10 



Examining the quantum properties of the cut-circle billiard in more detail, we show how five of the energy levels of 
the cut billiard behave as the cut size is varied from uj = 50° tow = 70° in Fig. 0. The states shown are continuations 
of the circle cavity degenerate pairs (I = 3, n = 2), (I = 6, n = 1), and one of the pair, (I = 1, n = 3). 

Here and in later numerical computations, we are using dimensionless variables. The dimensionless energy e is given 
by 

e =?£ E =MU£ (20) 

where n — kR is a dimensionless wavevector. To get an idea of how this translates to the physical dimensions of our sys- 
tem, consider an electron density of n«4x!0 11 cm~ 2 , which is a typical value measured for a two dimensional electron 
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gas in a GaAs/Al .3Gao.7As heterostructure |0. The corresponding Fermi wavevector is kf = y/2wnm 1.6x 10 8 m 1 , 
the Fermi energy Ej = ttH n/m w 0.014 eV (for an effective electron mass in GaAs of m — 0.067m e ). This yields for 
e = 100, which is the maximum energy we used in the computations, R rs 90 nm. 

In Fig. we observe a (distant) avoided crossing, which occurs between states 2b and 3a in the interval lu = 53° 
to oj = 63°. We also sec an actual crossing in this plot. This occurs between the states labeled 2a and 2b at about 
w = 65°. 

However, as discussed in Ref. fl2]| , separate eigenvalues of a generic Hamiltonian can in general only be brought to 
coincide if at least three parameters are varied. Thus, at first sight, we would expect another level repulsion instead 
of the actual crossing, since we only vary u>. Yet we ignored that the Hamiltonian of the cut-circle is not completely 
"generic" . Even if the cut is inserted, the system retains parity symmetry with respect to an axis perpendicular to 
the cut, i.e., 

P c V(j>,tp)=*(p,{-<p)=¥(p,<p), (21) 

where is the parity operator for a symmetry axis through the center of the cavity, at angle C (i.e., the angle at 
which the center of the cut is placed). Because of parity symmetry, we have two different symmetry classes, namely 
the classes formed by even and odd states. Since parity is a discrete symmetry, states of one class cannot interact 
with states of the other, thus both classes are completely independent of each other. Therefore, crossings between 
states of different parity may occur. Only within each parity class, i.e., among states with the same parity, crossings 
remain impossible. 

In Fig. |3|, we show energy eigenstates corresponding to the eigenvalues 2a, 2b, and 3a in Fig. ||, for three different 
cut sizes taken before the avoided crossing, between the avoided and actual crossings, and after the actual crossing. It 
can be seen that 2a has even parity, and 2b and 3a both have odd parities. This explains why 2a and 2b can undergo 
the crossing, whereas a crossing between 2b and 3a is avoided. 

We also see that the state 2a, which does not participate in the avoided crossing and undergoes the crossing with 
2b, remains unchanged. However, states 2b and 3a, which do undergo an avoided crossing, become mixed and loose 
their original identities. The mutual parity of states 2b and 3a allows them to couple, whereas this is impossible for 
states 2a and 2b with differing parities. 



IV. WAVEGUIDE WITH A CIRCLE CAVITY 



If we attach leads to the circle billiard, we have a matter waveguide and the eigenstates of the billiard become 
unstable and may decay. These unstable states strongly affect the scattering dynamics of the waveguide. They 
appear as resonances in the transmission probabilities through the waveguide, 

Z> = X;M 2 . (22) 

a 

for incoming mode (3 in lead 1, where a runs over the outgoing modes in lead 2. We can examine the individual 
Tps or consider their sum over all possible incoming modes in lead 1, which is proportional to the Landauer-Biittiker 
conductance p3[ , 

G=^£7>. (23) 



As shown in Ref. ||, there is an energy range where, for conical leads, mode (3 has nonzero transmission, whereas 
for straight leads with the same opening angles it would still be evanescent. This occurs for energies too small for 
the /3th transverse standing wave to fit into the lead opening at the junction to the cavity, such that the wave has 
to "tunnel" into the cavity. The condition for this tunneling (see Ref. ||) is k < (3ir/di, if rfj — 2i?sin(A#i/2) is the 
width of the lead opening (for opening angle A9i). Thus, tunneling occurs for energies such that 

. _ P 2 TT 2 R 2 _ f3 2 TT 2 

€<ep ~ 2d 2 ~ 8sin 2 (A^/2) ' ( } 



where e is the dimensionless energy defined in Eq. (|2C 

In addition to the transmission probabilities and Landauer-Biittiker conductance, it is useful to look at the Wigner 
delay times, 
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= , (25) 

where 9 n are the phases of the S-matrix eigenvalues s n — e l9n . The Wigner delay times are a measure of the time 
delay of the electron due to the presence of the cavity. In the following plots, we will use dimensionless delay times 

80 Tit w 

T " = IT = -%» ' ( 26 ) 

(For R ps 100 urn, r„ « 20, and to* = 0.067to 6 , the Wigner delay time is of the order of 0.1 ns.) From Ref. it is 
clear that peaks in the Wigner delay time correspond to resonances and poles of the scattering matrix in the complex 
energy plane. Therefore we can use peaks in the delay time spectra to indicate the resonances of our system. 



A. Symmetric placement of leads 

In Fig. ^, we show numerical data obtained for an electron waveguide with a circle cavity with symmetrically placed 
leads, one at 6\ = 0° and the other at 62 = 180°. Both leads subtend the same angle A(9i = A82 — 20°. The 
dimensionless energy e is varied in steps of 0.5. 

The transmission probability in Fig. 0(a) shows the quantities Tp for incoming modes (3 — 1,2 in lead 1. The 
first mode {(3 = 1) starts transmitting (i.e., becomes significantly larger than zero) at about e ss 20. The second 
mode starts transmitting too closely below 100 to be visible in the plot. The transmission for higher modes is zero 
everywhere in the energy interval below 100, and therefore these modes were not included in the computations or the 
plot. 

The threshold energies ep, for which the modes would start propagating in the case of straight leads with A9 = 20°, 
are ei = 40.91, and I2 — 163.66 for the first two channels. This means that for the conical leads, we observe tunneling 
of the first mode in the energy range between e sa 20 and 40. 

Examining the correspondence between eigenenergies of the closed circle (indicated by filled circles for degenerate 
levels and open circles for non-degenerate levels) and resonances of the open system, one finds that there is obviously 
a strong association. In energy intervals without eigenenergies, the transmission changes smoothly, whereas there are 
sharp changes close to eigenenergies. In the tunneling regime, we see peaks close to the eigenenergies. The peaks tend 
to be slightly shifted to the left. Above the tunneling threshold, we also observe valleys, e.g., at e = 67.5,74.7. In 
this regime, most peaks and valleys in the vicinity of an eigenenergy are too close together to tell which is associated 
with the eigenenergy. Since the lead openings themselves alter the geometry of the system, we cannot expect the 
resonances to be found exactly at the same positions as the eigenenergies. 

In Fig. ^(b), we plot the dimensionless delay times r„. The energy derivatives at each point are approximated by 
carrying out two computations for each point, with energy spacing Ae = 10~ 5 . Since the transmission probabilities 
Tp vanish for (3 > 3, we only have to consider the two channels (3=1 and (3 = 2. Thus, the S-matrix is a 4x4 matrix 
and we get four eigenphases and four delay time curves. 

Here, unlike for the transmission probabilities, we only need to look for peaks in the Wigner delay times at the 
resonance energies, in the tunneling as well as in the normal regime. The first two time delays display their first peak 
around e w 15, the others at e « 75. At these energies, the corresponding transmission probabilities are still too 
close to zero to produce a visible peak in Fig. ||(a) . This shows that the delay times are a more sensitive indicator of 
resonances than the Tp's. 

We see that the degenerate eigenenergies of the circle billiard do not produce double resonances when leads are 
attached at the special angles 9\ = 0° and 62 = 180°. Every eigenenergy of the circle billiard is associated only with 
a single delay time peak. Although placing leads to the circle cavity destroys its radial symmetry, the breaking of the 
degeneracy is not observed in this case. 



B. Asymmetric placement of leads 

In Fig. ^, we show the same waveguide, but with asymmetric placement of leads instead of the symmetric placement 
we had in Fig. ||. Lead 1 is placed at Q\ = 0° and lead 2 at 6 2 = 125°. 

Now we not only see again the close correspondence between delay time peaks and eigenenergies of the circle, but 
we also observe double peaks at most of the degenerate billiard eigenenergies, most evidently at e =24.6, 28.8, 35.4, 
38.5, and 49.4. At higher energies, the observation of double peaks is more difficult because the second mode starts 
transmitting. (Some of its first peaks are too small to be seen in the plot, but can be found by looking at the data 
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file.) Many of the double peaks show a small energy spacing: One peak is slightly shifted to smaller energies, the other 
one to higher energies. Here, the degeneracy is obviously broken by the addition of the leads. For the non-degenerate 
eigenenergies on the other hand, no splitting can be observed. Especially at e = 37.4, this becomes very obvious. 

The difference between this case and the previous case with symmetric lead placement (6*2 = 180°) is that, for 
symmetric placement, we had two discrete symmetries, namely invariance under the parity transformations Pqo and 
Pgo° [where the parity operators are defined as in Eq. (|2l|)l . When the leads are placed asymmetrically, the symmetry 
under Pqo, with respect to the axis bisecting both leads, is destroyed, and we are only left with one parity symmetry 
under Pg2 5° ■ The breaking of one parity symmetry seems to be the reason for the splitting of resonances seen in 
Fig. 1(b).' 



C. Numerical accuracy 



To get an estimate of the numerical error, the sum of the absolute values squared of all reflection and transmission 
amplitudes was calculated. Due to the unitarity of the S-matrix, this value has to be the dimension of the S-matrix, 
which is 4 here. The relative error, 

i - 1 E (M 2 + M 2 + l^l 2 + iCpl 2 ) , (27) 

a, (3=1 

was usually smaller than 0.2% for the individual data points. The errors of the absolute values of the S-matrix 
eigenvalues (which should be 1) were of the same order of magnitude. For the delay times, there is no easy way to 
estimate the error. 

The most crucial factor determining numerical accuracy was found to be the number 2n 7 of para mete rs (exp ansio n 
coefficients r a p,t a p, b v p) used to approximate the total wavefunction <!/ of our system [see Eqs. (|T^), ( A6), and ( A13 )]. 



An estimate for a suitable choice of this number can be found by comparing the Fermi wavelength, A/ — 2n/kf, to a 
typical length scale of the system, e.g., the length of the boundary ps 2itR (for zero or small cut size). If we assume 
that the total wavefunction of our system varies rather smoothly on a length scale of Ay, it seems plausible that 
a "small" number of parameters, say about 10, should be sufficient to approximate '5 on an interval of length A/. 
Therefore, we expect that such a "small" multiple of n — kfR = v / 2e="the number of A/-intervals needed to cover 
the cavity boundary" should give a reasonable order of magnitude estimate for 2n 7 . 

In Fig. ^|, we compare the transmission spectra [Fig. ^(a)] and the spectra of the time delay sums [Fig. ||(b)] for 
different values of n 7 . As deviations are expected to be most pronounced for large cut size, the calculations shown 



were done exemplarily for the symmetric cut-circle waveguide to be examined in Sec. VA, at ui — 80°. 

We find noticeable deviations for n 7 = 30 and, to a small extent, also for n 7 = 60. The values for n 7 = 90 
practically coincide with those for n 7 = 120. Therefore, in view of computation time, we used n 7 = 90 in all 
computations presented in this paper. 



V. WAVEGUIDE WITH A CUT-CIRCLE CAVITY 

We now consider a waveguide with a cut-circle cavity. Fig. [7] shows the sum of the partial delay times, t ~ 5Z n=1 t„, 
and the Landauer-Biittiker conductance G. For this computation, leads and cut are placed asymmetrically. The leads 
are centered at B\ = 0° and 82 = 166°, and the cut at C = 225°. The cut size u> varies from 0° (circle) to 80°, with 
step size 5°. The dimensionless energy e is varied in steps of 0.25. 

In the delay time spectrum (upper plot), we observe chains of separated double peaks, corresponding to the split-up 
branches of degenerate circle eigenenergies. If we compare the positions of maxima in the Wigner delay time to the 
positions of eigenenergies of the closed system, we find a good agreement between (open system) resonances and 
(closed cavity) eigenenergies. 

The conductance (lower plot) shows similar features. However, the first visible chains of resonances occur at 
significantly higher energies than in the delay time plot. For energies between 20 and 40, both plots show basically 
the same structures. Above the first tunneling threshold at t\ = 40.91, the conductance reflects the onset of "normal" 
(compared to tunneling) transmission of the first mode. In this energy regime, we cannot tell any more whether 
resonances are indicated by peaks or dips in the conductance. 

We will now compare the properties of the cut-circle cavity waveguide for symmetric and asymmetric lead place- 
ments. Note, however, that although this may seem to be analogous to the analysis from the previous section (for the 
full-circle cavity), the insertion of the cut already leaves us with only one parity symmetry (under P^). As a result, the 
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"symmetric" cut-circle waveguide is analogous to the asymmetric full-circle waveguide rather than to the symmetric 
full-circle waveguide. Finally, in the "asymmetric" cut-circle waveguide, there are no geometric symmetries present 
at all. Since this could not be achieved by merely placing two leads to a circle, this case has no parallel in Sec. [V. 



A. Symmetric placement of leads and cut 

Let us examine the Wigner delay times for the energy range where we found an actual crossing and an avoided 
crossing in Fig. ^ for the closed cut-circle billiard. We first consider a system in which the leads are placed at 9\ = 0° 
and 6*2 = 120°, symmetrically with respect to the cut at £ = 240°. The lead openings are A6>i = A62 = 16°. The plot 
of the sum r of the Wigner delay times is shown in Fig. ^. 

We observe five delay time peaks, which are Fano resonances resulting from the closed cut billiard eigenstates la, 
lb, 2a, 2b, and 3a shown in Fig. ^. Starting from the left, we label these delay time peaks in the same manner. For 
cut size u = 20°, we designate the delay time peaks, going from left to right, as lb, 2a, 2b, 3a, and 3b. la is the 
resonance chain in the upper left corner. The first two delay time peaks, la and lb, are Fano resonances resulting 
from the circle eigenstate with (I = 3, n = 2), the following two from (I = 6, n — 1), and the last two, 3a and 3b, from 
(Z = l,n = 3). 

Again, we observe an avoided crossing between 2b and 3a, and the two delay time peaks 2a and 2b actually cross. 
Like for the closed cavity, the crossing is possible due to the parity symmetry of the system, which has been preserved 
by adding the leads symmetrically with respect to the cut. Since 2a and 2b have different parities, these eigenvalue 
chains are allowed to cross. 

The symmetric placement of leads has not broken any additional symmetries. Only the energies and cut sizes at 
which the avoided crossing and actual crossing take place are slightly shifted from the case of the closed cut billiard. 

B. Asymmetric placement of leads and cut 

In Fig. H we show the same Fano resonances as in Fig. ||, but for asymmetric placement of leads. Now we place the 
leads at 6>i = 0° and 2 = 166° and the cut at C = 225°. 

The crossing of resonances 2a and 2b, which occurred for the symmetric leads and for the closed cut-circle billiard, 
has now become an avoided crossing. (In fact we now have a sequence of two closely spaced avoided crossings, similar 
to the case studied in Ref. juj.) Here, the parity symmetry of the closed system is destroyed by placing the leads 
asymmetrically. Therefore, in the asymmetric open system, we no longer have two separate classes of states (namely 
states of even and odd parity) and the crossing, which was possible in the closed system and in the symmetric open 
system, is avoided here. 

C. Cavity wavefunctions 

We now look at the distribution of electron probability in the cavity at the resonance energies. This is shown in 
Fig. [l^ for the asymmetric system. A three-way avoided crossing has mixed the resonance states, and has destroyed 
the identity they had for smaller cut size. The upper row displays the four states before the avoided crossing, which 
still reflect the quantum numbers of the associated circle states very well. For the middle row, the character of states 
2a, 2b, and 3a has been changed. Only the resonance state lb, which never comes really close to the crossing region, 
has preserved its characteristic structure (two rings with six peaks each) quite well. In the bottom row, after the 
second avoided crossing, 2a and 2b have mixed even further. 

VI. CONCLUSION 

In this paper, we investigated a waveguide with a cut-circle cavity by calculating conductances and Wigner delay 
times. We attached conic leads to the cavity in two ways, symmetrically and asymmetrically. We observed two kinds 
of avoided crossings when the leads were placed asymmetrically, breaking the parity symmetry of the closed cavity. 
One avoided crossing was due to the chaos of the closed cavity, i.e., due to level repulsions between states in the 
same symmetry class, which are known to occur for fully chaotic systems. The other avoided crossingwas due to the 
breaking of the discrete symmetry (the parity in our case) by placing leads asymmetrically (see Fig. g). 
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There is now great interest in exploring the effects of underlying chaos on the dynamics of open quantum systems. 
However, open quantum systems must make contact with the outside world. This contact itself may induce some of 
the signatures of chaos. This fact has also been studied extensively by Jung and Seligman Q for classical scattering 
systems. The avoided crossing between 2a and 2b, shown in Fig. provides another example of this effect. 



APPENDIX A: NUMERICAL METHOD 



1. Basic concept 

In order to obtain an explicit expression for the S-matrix, we need to find the stationary solutions 'J of Eq. (p]). 
This can only be done numerically for our system. For that purpose, we use the boundary integral method described 
m Ref. @. Below, we will give an outline of this method adapted to our particular problem. 

The boundary integral method is based on the use of Green's Theorem, which, for two functions / and g, states 
that 



(Al) 



if r is a closed contour confining an area A, I is the longitudinal coordinate along T, and n is the normal coordinate. 
For the boundary integral method, we choose for the functions / and g the wavefunction ^ and a weight function <!> 
that satisfies the Hclmholtz equation 

V 2 $ = -fc 2 $ . (A2) 

Then, since Vfr i s a solution of the Schrodinger equation [Eq. (|l|)], i.e., \& also satisfies V 2 ^ = — fc 2 \I>, the right hand 
side of Eq. ( |A1[ ) is identically zero. Therefore we get 

*( , 5T-5r , )-° < A3 > 

as the basic equation of this numerical approach. (Note that the difference between ^ and $ is that there are 
no boundary conditions for <f>, whereas \& has to satisfy the boundary conditions imposed by the geometry of the 
waveguide.) 

For our problem, we choose T to follow the walls of the cavity. For the places where the leads come in, T is taken 



to be the continuation of the circular wall segments in order to construct a closed contour. Then, Eq. (A3) allows 
us to obtain the S-matrix and the probability distribution inside the cavity for any incoming particle energy. This is 
done by inserting (approximate) expansions of W with, say, n cuto g yet u nkn own expansion coefficients. As we need 



^cutoff equations to solve for these expansion coefficients, we apply Eq. (A3) with sufficiently many different weight 
functions <&„ (n = 1, . . . , n cuto ff)- 

In principle, there are several choices possible for the $„'s, as long as they satisfy Eq. (|A^). For numerical 
convenience, we used the same functions as for the expansion of \& inside the cavity, namely 4> 7 = / 7 with the / 7 's 
given by Eq. (||), where 7 = — n 1 + 1, . . . , n 7 (i.e., n C utoff = 2n 7 ). For n 7 , we also used the same value as in Sec. [I A, 
n 7 = 90, which allows us to determine n cuto ff = 180 expansion coefficients. 



2. S-matrix 



In order to compute the S-matrix, we split up Y into parts C\,C% across the lead openings, PxyP^P^ along the 
circular walls, and Pq along the cut (see Fig. pi] ). 

We first regard a part Iq of the integral in Eq. (A3), for which we only integrate over the lead openings C\ and 
C2 instead of T. Inserting for the / 7 's and for \1/ its expansion into lead channels from Eq. (|l^) (for an electron 
incident in channel (3 in lead 1), we obtain 



Ic 




dn 



f- 
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Oo )l0 - ( Z L hf) + £ Ta P [( Z + )ia ~ ( z L + )l 
a 

+ y j tap [(zq + ) 7Q — (2:^)70] , 



(A4) 



where we introduced the abbreviations 



dn h ' 



(' - " i - : / dl Y i=t ^ 



(A5) 



with i =1, 2 denoting lead 1 and 2. (For an electron incident in lead 2, the expressions look similar.) 

Along the walls of the cavity, the electron wavefunction is zero, but it may have a finite slope perpendicular to the 
wall. We can expand this slope in terms of a complete orthonormal basis along the wall. Therefore, we can write 



*/3,Pi = : 



dn 



(A6) 



with some set of basis functions in variable I. The index i = is used to denote the cut, and i = 1,2,3 denote the 
three different circular wall segments. We are free to choose any complete basis, and for the circular wall segments 
we choose a Fourier basis, 



e^(0 = A/^cos 



AT 
2 , 



(A7) 



where 8i is the angular position of the center of wall segment i = 1,2, 3, Aj is its opening angle, v = 0, 1, . . ., and 
St — < Z < 5j + For the cut, we could use the Fourier basis, too. However, as this results in very long- 
computation times, a basis of small triangles proved to be more efficient. Thus, along the straight wall we use 

f+1/2 



6 0) d) 



-TT for — ^TT < I < ^44 



i Z - J for ££±£ Z < Z < Z 
else 



(A8) 



for < Z < Zo, where Zo = 2i?sin(o;/2) is the length of the cut, w is its opening angle, and v = 0, 1, . . . , n„ . Now we 
can write down the part Ip of the integral in Eq. ( A3), for which we only integrate over the wall segments Po, . . . 
instead of T. Inserting the expansions from Eq. ( A6), we get 



££«/ <^ 3 /, = 



(A9) 



i=0 1/ 



i=0 f 



with the abbreviation 



7^ 



dZ • 



(A10) 



p, 



Now we can combine these results to obtain a matrix equation f or t he unknown expansion coefficients, which can 
easil y be solved numerically. With the integrals Ic and Ip, Eq. (AS) becomes Ic + Ip = 0. Inserting Eqs. ( A4) 
and (A9) yields (after a slight reordering of terms) 



EE^V^ = - [(4 hf> ~ (4 hf>_ 



(All) 



or, in compact matrix form: 
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Z+-Z+, -X) t =(z L --z -) /3 . (A12) 



The matrices Zj, Zj, and X, as well as the column matrix (z L — z ) „ can be computed numerically from Eqs. ( |A5| ) 

and ( |A10| ). Note that in order to include both leads (i = 1,2) in Zq~ and Zj, the first columns of these matrices 
correspond to lead 1, and the columns for lead 2 are simply appended. (Similarly, X is composed of the different wall 
parts, i = 0, . . . , 3, also columnwise.) 

The quantity I t I is a column matrix containing the unknown transmission and reflection coefficients for a 



b 







(i) 

particle incident in channel p in lead 1, as well as the coefficients b v p, which describe the normal derivative of ^fp 
on the walls. As this is the quantity which we want to solve for, we need to invert the matrix ( Z q — Z \ , —X ) . 
Therefore, this matrix has to be square. This can be achieved by using appropria te c utoff numbers n ai (i = 1,2) in 



the lead expansions [Eq. (16)] and n Vi (i = 0, . . . , 3) in the wall expansions [Eq. (A6)]. Since the matrix consists of 
2n 7 rows, these numbers have to satisfy the relation 

2 3 

Ua * + J2( n ^ +1 ) = 2n i ■ ( A13 ) 

i=l i=0 

(Note that v = 0, . . . , n Vi , therefore we actually sum over n Vi + 1 terms for the wall parts.) For best numerical accuracy, 
the relations of the cutoff numbers are chosen approximately equal to the length relations of the corresponding 
integration paths. 



3. Cavity probability amplitude 

Once the elements of the S-matrix have been obtained by the method described above, it is an easy task to calculate 
the probability amplitude for the electron state inside the cavity. We know that the wavefunction, which, inside the 
cavity, is described by Eq. (|J), is zero along all the walls and that is continuous at the interface between the leads 
and the cavity. Thus we require that along the walls, 

]Tc 7/3 />,^) =0, (A14) 

7 

and along the interface C± and C% between the leads and cavity 

^c 7/3 / 7 (p, <p) = 9 p (p, ip) , (A15) 

7 

where (in the leads, for incoming mode (3) is known from the previous section. It is straightforward to solve for 
the c 7) 3's, e.g., by evaluating the above expressions at 2n 7 different points on T, and inverting the resulting system of 
linear equations. 
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FIG. 1. Geometry of the scattering system with cavity radius R, lead positions at 6i and 82, lead opening angles A9i and 
A02, cut position at £ and cut angle lo. 
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FIG. 2. Five energy eigenvalues of the closed cut billiard as a function of dimensionless energy e and cut size lo. From left 
to right we label them as la, lb, 2a, 2b, and 3a. For lo — 0°, the states la and lb are degenerate with quantum numbers 
(I = 3, n = 2), the states 2a and 2b are degenerate with (I = 6, n = 1), and 3a is one of a degenerate pair with (I = 1, n = 3). 
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FIG. 3. Three energy eigenstates corresponding to energy eigenvalues 2a, 2b, and 3a in Fig. g. The cut is located at £ = 270°. 
The top row has lu = 50° and consists of (a) State 2a for e = 50.2, (b) State 2b for e = 50.5, and (c) State 3a for e = 51.9. The 
middle row has u) = 61.4° and consists of (d) State 2a for e — 51.0, (e) State 2b for e = 51.3, and (f) State 3a for e = 52.4. The 
bottom row has ui = 68° and consists of (g) State 2b for e = 51.5 and (h) State 2a for e = 51.7. Note that, from left to right, 
we show states with increasing energy. This means that 2a and 2b switch position from the middle to the bottom row, as their 
eigenenergies cross. 
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FIG. 4. (a) Transmission probability and (b) delay time spectra for lead openings A9\ — 20°, A#2 = 20° and lead positions 
6>i = 0° and 82 = 180° . Eigenvalues of the closed circle are marked by filled (degenerate eigenvalues) and open (non-degenerate 
eigenvalues) circles. 
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(a) t„ 




FIG. 5. (a) Transmission probability and (b) delay time spectra for lead openings A9i — 20°, A62 = 20° and lead positions 
81 — 0° and 62 = 125°. Eigenvalues of the closed circle are marked by filled (degenerate eigenvalues) and open (non-degenerate 
eigenvalues) circles. 
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FIG. 6. (a) Transmission spectra and (b) spectra of the sum of partial delay times in the dimensionless energy interval 
47.5<e<52.5, for cut size u><80° . The lead openings are A61 = A02 = 16°. Lead and cut positions are 9\ — 0°, 82 = 120°, and 
£ = 240°. The results for different values of n 7 are compared. 



15 




FIG. 7. (a) Sum of the partial Wigner delay times and (b) conductance for lead opening Ad = 20° and varying cut size ui. 
The leads are centered at 9\ = 0° and 6 2 = 166°, and the cut is centered at £ = 225°. 




FIG. 8. Spectrum of the sum of partial delay times for the energy interval 47.5<e<52.5 and cut size interval 20 o <cj<80°. 
The lead openings are Adi = A82 = 16°. Lead and cut positions are 9i = 0°, 62 = 120°, and £ = 240°. 
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FIG. 9. Spectrum of the sum of partial delay times for the energy interval 47.5<e<52.5 and cut size interval 52°<lj<80°. 
The lead openings are A#i = A62 = 16°. Lead and cut positions are 9i — 0°, 62 = 166°, and C, = 225°. 




FIG. 10. Electron probability distribution in the waveguide cavity at the resonance energy for Fano resonances shown in 
Fig. ^[ The lead openings are A#i = A62 = 16°. Lead and cut positions are 0i = 0°, 62 = 166°, and £ = 225°. The top row 
has lu — 50° and consists of (a) State lb for e = 47.8, (b) State 2a for e = 49.2, (c) State 2b for e — 50.2, and (d) State 3a 
for e = 51.1. The middle row has uj = 68° and consists of (e) State lb for e = 48.5, (f) State 2a for e = 50.7, (g) State 2b for 
e = 51.3, and (h) State 3a for e = 52.3. The bottom row has u> = 78° and consists of (i) State lb for e = 49.2, (j) State 2a for 
e = 50.7, (k) State 2b for e = 52.2. 
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FIG. 11. Integration contour for the cut-circle waveguide, consisting of Ci, Ci, Pq, Pi, Pi, and P3. 
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